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ABSTRACT 

Numerical simulations are performed on different lattice sizes of chiral U(l) 
and SU(2) scalar-fermion models with explicit mirror pairs of fermions in the 
broken symmetry phase. Relevance of these models to the electroweak theory is 
discussed. Shift symmetry of the action is exploited to decouple the unwanted 
mirror-fcrmion. Upper and lower bounds on the renormalized scalar mass as a 
function of the renormalized fermion mass are obtained in the SU(2) version of 
the model. Our numerical data are found to be in qualitative agreement with 
one-loop results. No evidence for a nontrivial fixed point has been observed. 
The mirror- fermion models are likely to be trivial in the continuum limit. Our 
Monte Carlo data show that with a 180 GeV top quark, the Standard Model 
Higgs boson should have a mass between 100 and 800 GeV. 

1. Introduction 

It is believed that the only possible nonperturbative (i.e. strongly interacting) 
sectors in the electroweak theory of the minimal Standard Model are the Higgs and 
heavy fermion sectors. (Here heavy or light is meant to be compared to 100 GeV scale, 
the energy scale of the weak interaction.) If the electroweak theory is trivial in the 
infinite cutoff limit, then we can obtain upper and lower bounds on the Higgs mass 
as a function of the heavy fermion mass. Since those sectors can become strongly 
interacting, it is important that nonperturbative investigation be carried out. In this 
report, the nonperturbative method is Monte Carlo simulation on the lattice. 

Because there are lots of bare parameters in the electroweak theory, we need to 
make some approximations before lattice calculations become feasible. The approx- 
imation we make is to ignore all the gauge couplings assuming that all gauge cou- 
plings produce perturbative effects. (Although QED has a Landau pole, its coupling 
at 100 GeV scale can still be safely ignored because Landau pole presumably appears 
at some astronomically high scale.) We call this the zeroth-order approximation. We 
also ignore all the light fermions. Therefore, the only fermion left is the top quark 
(and/or the fourth family fermions should they exist). 

It is well known that naive fermion action on the lattice suffers from species 
doubling. The ways to handle fermion doublers are the Wilson fermion and the 
staggered fermion a. In this report, we take the approach of the Wilson fermion. 
However, the electroweak theory is a chiral gauge theory. Even after we ignore the 
gauge couplings, we still have global chiral symmetries left. The naive Wilson fermion 



action breaks the chiral symmetry explicitly and is not good for the electroweak 
theory. In order to cure this, two models were proposed. One is the Smit-Swift 
model i. The other is the mirror-fermion model i. We will report on our lattice 
study of the electroweak theory using the mirror-fermion model in the zeroth-order 
approximation. 

The global symmetry of the mirror-fermion model can be U(l) or SU(2). It is 
found that the U(l) and SU(2) versions have the same qualitative properties. Since 
SU(2) symmetry is closer to the electroweak part of the real world, we will concentrate 
on the SU(2) mirror-fermion model. 



2. Lattice Action and Decoupling Limit 

The lattice action of the mirror-fermion model is a sum of the 0(4) (= SU(2)l 
SU(2) R ) symmetric pure scalar part S v and fermionic part S^: 
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(f x is the 2 ® 2 matrix scalar field, and *f? x = {ip x ,Xx) stands for the mirror pair of 
fermion doublet fields (usually ip is the fermion doublet and \ the mirror fermion 
doublet). In the usual normalization conventions for numerical simulations we have 
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Here the two fermion fields ip and \ are Wilson fermions, K is the fermion hopping 
parameter, r the Wilson-parameter, which will be fixed to r = 1 in the numerical 
simulations, and the indices L, R denote, as usual, the chiral components of fermion 
fields. One can easily see that the above action has an explicit chiral symmetry in 
the sense that left- and right-handed components of the fermion fields can transform 
independently. However, the right-handed x-held has to transform exactly as the 
left-handed ^-field, and left-handed x-component as the right-handed ^-component. 
(So an off-diagonal bare fermion mass term can exist in the action.) Therefore, ip and 
X are mirror partners of each other. In our simulations, we take normalization that 
fermion-mirror-fermion mixing mass is fi^ x = 1 — 8rK. 



The fermionic part is given here for a single mirror pair of fermions. Taking the 
adjoint transforms fermions to mirror fermions and vice versa, but as noted before, 
without SU(3) co iour®U(l)hy P ercharge gauge couplings they are equivalent to each other. 

All renormalized quantities can be defined. (See @ for details.) Since the symmetry 
of the action does not exclude the mixing fermion mass term, we will have a nonzero 
renormalized fermion mixing mass [Ar. The physical fermion states will be linear 
combinations of the renormalized ip and \ fields in general. One can always tune the 
bare parameters to have a zero renormalized mixing. It turns out that in most cases 
this tuning is rather difficult and time consuming in the simulation. 

So far we have not discovered mirror-fermions in the experiments, we do not want 
to have them in the physical spectrum. There are different ways to decouple those 
unwanted fermions. We can tune the bare parameters such that mirror-fermion has 
its mass of the order of the cutoff scale. Hopefully this heavy mirror-fermion will 
decouple. The more elegant way is to take the advantage of the symmetry of the 
model and is explained below. 

The action shown in eq.(|2|) is invariant under the following transformations of the 
fermion fields 

ipx — > i/> x + u , 4>x — » i> x + v if Gip = ; (3) 

Xx -> Xx + e , Xx -> Xx + e if G x = (4) 

where u, e do not depend on x. Due to the above symmetry (called shift symmetry 
from now on), one can easily derive the corresponding Ward identity and show that 
Gr x = identically at G x = 0. Furthermore, from the Ward identity for the fermion 
propagator, one can show that [ir will vanish identically if we choose K — 1/8 in 
addition to G x = 0. This analytic result tells us that at G x = 0, K = 1/8, there 
is no mixing between the fermion and mirror-fermion, and the renormalized Yukawa 
coupling of the mirror-fermion is zero. Therefore, the mirror-fermion is massless and 
is not coupled to the scalar and fermion fields. In other words, at G x = and 
K = 1/8, the mirror-fermion behaves like the massless right-handed neutrino and is 
decoupled from the physical world i. We call the above choice of the bare parameters 
the mirror-fermion decoupling limit. 

At the energy scale of the weak interaction (around 100 GeV) , all gauge couplings 
are weak. We therefore take the zeroth-order approximation by ignoring all gauge 
fields. We think that the only possible sources of nonperturbative physics in the 
electroweak theory are the Higgs and heavy fermion sectors. Notice that the shift 
symmetry of our mirror-fermion model is destroyed if the gauge fields are put in. 
Mirror-fermion will be coupled to the physical world via gauge interactions even if 
we set G x = 0, K = 1/8. We simply assume that the zeroth-order approximation 
reported here serves as a good approximation to the electroweak theory without 



worrying how to decouple the mirror-fermion in the full theory. 



3. Vacuum Stability and the /3-Functions 

The upper bound on the Higgs mass comes from the notion of triviality §. As for 
the lower bound, notions of triviality and vacuum stability are both required i. 

In order that a quantum field theory is self-consistent, its energy spectrum has 
to be bounded from below. Or we say that its vacuum has to be stable. Since the 
minimum of the effective potential corresponds to the ground state (i.e. vacuum) 
of the theory, a self-consistent quantum field system should have an effective poten- 
tial which is convex. Namely, as the field strength increases, the effective potential 
cannot go down to negative values indefinitely. In terms of the language of renormal- 
ization group, the running scalar coupling constant cannot go to negative values as 
we increase the energy scale @. 

What describe how the couplings change as the energy scale changes are the (3- 
functions. The exact /3-functions and effective potential are usually unknown. People 
rely on loop expansions to derive them. However, loop expansions can break down 
when the system becomes strongly interacting, we must go beyond perturbation the- 
ory. At the nonperturbative level, we find that the partition function of the Euclidean 
mirror-fermion model simply blows up when A is negative. If this is so, then the ef- 
fective potential is not defined because the theory is not self-consistent. In order that 
the theory exists, A has to be positive. Then the partition function exists and the 
effective potential can be defined. Since the exact effective potential is defined via a 
Legendre transformation, its shape will be convex and has a minimum once it exists. 
Hence, at the nonperturbative level, once A is positive, the theory will have a stable 
vacuum. 

If the theory is trivial in the continuum limit, then its exact /^-functions will 
behave qualitatively like the one-loop ones. There we can see that at some large 
values, /3-function of g can turn negative. This shows that g can become negative as 
we go to higher energy scales, thus spoiling the stability of the vacuum. In order that 
the vacuum state be stable, gn cannot be too small. Hence, we can obtain a vacuum 
stability lower bound on gp as a function of Gr^. Since the smallest value A can take 
without destroying vacuum stability is zero (or + ), the lower bound on gp naturally 
corresponds to a very small but positive A @. 

Thus in lattice study of the SU(2) mirror-fermion model, we can first decouple the 
unwanted mirror-fermion by setting G x — 0, K = 1/8. If the mirror-fermion model is 
trivial in the continuum limit, Monte Carlo data on gn and Gp^ obtained at A = oo 
and A = + will give us upper and lower bounds on the Higgs mass as a function of 
the heavy fermion mass respectively. It goes without saying that both bounds are 
cutoff dependent. 

We also carry out one-loop renormalized perturbation theory on the lattice and 



Fig. 1. Higgs mass bounds as a function of heavy fermion mass. Curves are from one-loop pertur- 
bative calculations in the continuum at £ = 1.0 and 1.2. 

calculate /3-functions accordingly. We therefore can estimate the finite size effects and 
finite cutoff effects. 

4. The Simulation and Data 

The phase structure of the model was first explored. (See i for details.) We 
find that the phase boundary between the symmetric (PM) phase and the broken 
symmetry (FM) phase is a second-order phase transition line where the cutoff can be 
removed. 

Numerical simulations were performed on 6 3 ■ 12 and 8 3 • 16 lattices at A = 1CT 6 
and oo. The small positive value of A was chosen to mimic A = + . The Yukawa 
coupling G x was kept at zero, the fermion hopping parameter K was kept near 1/8 
on the finite lattice for exact decoupling of the mirror doublets !'!. Since it is the 
FM phase that is physically relevant, k and G$ are tuned such that we approach 
criticality from within the FM phase and have a reasonable cutoff to reduce finite size 
effects. Since Ward identities give us some analytic results for renormalized quantities 
at G x = 0, we put those results into our Fortran program to "guide" our Monte Carlo 
simulations. This helps to reduce the fluctuations. 



5. Conclusions 



Our numerical data are presented in Fig. 1. At all points of the simulation, we 
find that all fermion doublers are decoupled as expected. 

It is seen that our Monte Carlo data are in qualitative agreement with one-loop 
perturbative results. This may not be so surprising since tree unitarity (which requires 
probability conservation) upper limit on G\^ is 2tt. Below this value, the system 
should be weakly interacting. So far, no evidence for an ultra-violet stable fixed 
point has been found. Since all qualitative features of the one-loop /^-functions are 
supported, our mirror-fermion models are likely to be trivial in the infinite cutoff 
limit. Assuming that the SU(2) mirror-fermion model is trivial in the continuum 
limit, we obtain, in zeroth-order approximation, upper and lower bounds on the 
Standard Model Higgs mass as a function of the heavy fermion mass. If the heavy 
fermion in our model is top quark with its mass around 180 GeV, then we claim at 
£ = 1.0 that Standard Model Higgs has its mass between 100 and 800 GeV. 

Although the observed qualitative behaviour is certainly consistent with the one- 
loop perturbative scenario implying the triviality of the continuum limit, one has to 
keep in mind that present simulations are done at relatively low cut-offs. In particular, 
the evolution of the couplings towards smaller values at decreasing cut-offs should be 
investigated in the future in order to make sure that the model is indeed trivial. 
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